Imaging ac losses in superconducting films via scanning Hall probe microscopy 
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Various local probes have been applied to understanding current flow through superconducting 
films, which are often surprisingly inhomogeneous. Here we show that magnetic imaging allows 
quantitative reconstruction of both current density, J, and electric field, E, resolved in time and 
space, in a film carrying subcritical ac current. Current reconstruction entails inversion of the Biot- 
Savart law, while electric fields are reconstructed using Faraday's law. We describe the corresponding 
numerical procedures, largely adapting existing work to the case of a strip carrying ac current, but 
including new methods of obtaining the complete electric field from the inductive portion determined 
by Faraday's law. We also delineate the physical requirements behind the mathematical transfor- 
mations. We then apply the procedures to images of a strip of YBa2Cu307_5 carrying an ac current 
at 400 Hz. Our scanning Hall probe microscope produces a time-series of magnetic images of the 
strip with 1 fj,m spatial resolution and 25 /is time resolution. Combining the reconstructed J and 
E, we obtain a complete characterization including local critical current density, E-J curves, and 
power losses. This analysis has a range of applications from fundamental studies of vortex dynamics 
to practical coated conductor development. 

PACS numbers: 74.78.Bz, 74.25. Sv, 74.25. Qt, 07.79.-v 



I. INTRODUCTION 
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Developing high-temperature superconducting films to 
carry high current densities with low dissipation has 
proved to be an enormous challenge, resulting in compos- 
ite materials that are inhomogeneous down to the atomic 
scalcfi For characterizing inhomogeneous current flow, 
magnetic imaging has emerged as an important tool. By 
measuring the magnetic field in a plane above a film's sur- 
face, one can invert the Biot-Savart law to reconstruct a 
two-dimensional map of current density, J, in the film, 
as described in Section IIV Al ^'^'"^ Such a measurement 
can be obtained by various methods including magneto- 
optical imagingf^i^ii, scanning Hall probe microscopjsSi^, 
and scanning superconducting quantum interference de- 
vice (SQUID) microscopy. ^"'^^ 

Other imaging techniques, such as 
potentiometry^ scanning laser microscopyi 
scanning electron microscopy^** can map out the elec- 
tric field, E, that arises from vortex movement or 
other changes in local supercurrent density. In these 
techniques, to generate measurable electric fields, the 
superconductor is biased with a current / slightly greater 
than the dc critical current Ic at which vortices start to 
flow. In principle, such measurements could be combined 
with magnetic imaging of the same sample, and together, 
J and E would provide a complete, spatially-resolved 
electrical characterization of the material, including the 
local critical current density Jc, and the local power 
input, which can be calculated as J • E. 

In this work, we demonstrate that for I < Ic, time- 
resolved magnetic imaging can simultaneously determine 
both J and E in a superconducting film. Figure [T] sum- 
marizes the operations involved. The instantaneous mag- 
netic field B determines J, while the time rate of change 
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This relation only constrains the inductive portion of 
electric field, E^. To reconstruct the remaining electro- 
static part, Ep, we must impose the additional restric- 
tions that I < Ic, as discussed in Section |Vl and that 
E is parallel to J, discussed in Section IV CI However, 
many important applications that are not accessible to 
techniques operating above Ic do lie within these restric- 
tions, such as ac losses in superconducting films. 

Faraday's law has been applied previously to derive 
ac loss from a critical state model of magnetic fields in 
a homogeneous, infinitely long superconducting wirej ^^d^ 
In the present work, this method is reformulated to allow 
for the inhomogeneity of a real conductor revealed by 
our magnetic images. The method also applies to the 
case of magnetization decay due to fiux creep, which was 
recently analyzed using magneto-optic images.'^^ 

Here, we use a cryogenic scanning Hall probe micro- 
scope (SHPM) to obtain a series of images of B^, the 
component of magnetic field perpendicular to the sam- 
ple surface, as it evolves with time. The microscope is 
described in Ref. and further information is given in 
Section IITXI 

To analyze ac losses, we image a strip of the high- 
temperature superconductor YBa2Cu307_5 (YBCO), 
described in Section IIIBI We apply an ac current to 
the strip at 400 Hz, a typical operating frequency for 
applications,^'* and must image the magnetic response 
faster than this. Our scan speed — a few pixels per 
second — is slow in comparison, so instead of acquiring 
an entire sequence of images within one cycle of applied 
current, we obtain an average response over many cycles, 
as explained in Section III CI with 25 time resolution. 
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FIG. 1: A flow chart illustrating how we use a series of mag- 
netic fleld images to reconstruct current density and electric 
field in the sample, and finally power input. 



size leads to averaging of the field over 500 nm. Second, 
though its tip remains in contact with the sample during 
scanning, the sensor's slight tilt lifts the sensitive area 
1 /im above the sample surface. These factors yield 1 /im 
spatial resolution. 

The sensitive area is coated with a grounded gold 
gate to screen electric fields, followed by an insulating 
aluminum oxide layer to isolate the gate from sample 
voltagesii^ The oxide is also intended to provide protec- 
tion against mechanical wear. 

B. Sample: YBCO strip 

The YBCO film studied is 180 nm thick, grown epi- 
taxially by pulsed laser deposition (PLD) on a SrTiOa 
[001] substrate. Photohthographic patterning followed 
by argon ion milling removes parts of the film, leaving 
a bridge as shown in Fig. [2ja). The substrate is held in 
vacuum, attached to the microscope's copper coldfinger 
by a thin layer of low-temperature varnish. The coldfin- 
ger is cooled by a continuous flow of liquid helium. The 
film's Tc (defined by the maximum in dR/dT) is mea- 
sured in this cryostat to be 90 K. For imaging, the film's 
coldfinger is held at 40 K while the Hall sensor is held at 
54 K. However, the current return lead (to the right of 
the segment imaged) is narrower than the bridge (32 /im 
versus 50 /^m), and a magnetic scan (not shown) indi- 
cates that the applied current of 0.75 Arms exceeds its 
critical current. The consequent dissipation in the return 
lead may heat the sample several degrees above 40 K. 



Section IIIII presents this time series of magnetic im- 
ages, which are then transformed into images of current 
and electric field in Sections IIVI and [V] These quantities 
are combined in Sections IVII and IVIII to yield maps of 
dissipation and superconducting characteristics. 



II. EXPERIMENTAL PROCEDURE 
A. Apparatus: scanning Hall probe microscope 

The apparatus is described in Ref. H, but essential ex- 
perimental details and notable modifications are included 
here. The instrument generates magnetic images by ras- 
tering a Hall sensor above the sample surface, measuring 
magnetic field at each pixel. The sensor's Hall resistance 
is approximately proportional to the component of 
magnetic field perpendicular to the sample surface, with 
a field resolution of 4 x 10^^ G/VHz. 

The scanning stage, based on stepper motors driving 
micrometers, offers a 1 x 4 cm scan area for aligning to 
samples and macroscopic features, and can zoom in for 
200 nm positioning resolution. However, the image res- 
olution is limited by the sensor: First, its lithographic 



C. Data acquisition 

The images shown in Fig. [2Kc)-(j) are all derived from 
one scan as described in Ref. A 402.7 Hz ac cur- 
rent is applied continuously to the sample. The sensor 
rasters over a 100 by 150 /xm area with 0.5 /im pixel spac- 
ing. It pauses at each pixel and records the waveforms 
I{t) and B{t) for approximately 80 cycles of applied cur- 
rent. These cycles are overlaid and averaged. Values 
from the averaged waveforms are collated into images by 
their phase within the cycle. The complete set of images 
is presented as the movie BJE.avi in the supplemental 
materiaL-" The images correspond to 100 time slices over 
the cycle of applied current, each representing the aver- 
age magnetic field over a 25 /iS interval. 

III. MAGNETIC IMAGE RESULTS 

This work is primarily concerned with deriving quanti- 
ties of interest from time-resolved magnetic images, but it 
is worth first examining the images directly. From them, 
we extract a great deal of qualitative information that 
shapes the assumptions under which further quantities 
are derived. 
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FIG. 2: (Color online) Imaging current-induced flux penetration into a YBCO film: (a) Illustration of the sample geometry. 
The approximate area of the magnetic images is outlined by the brown box. Current is injected at the left and extracted 
from the downward-facing lead; the right-facing current lead and smaller voltage leads are not used (floating), (b) Applied 
current during a 0.75 Arms, 400 Hz cycle. Select times are marked, and the corresponding magnetic images are shown in 
(c)-(j). Image (c) indicates the scale and boundary of the experimental data (dashed green box); the background outside of 
the data is filled by a fit to a critical state model. The vertical dashed yellow line indicates the location of the cross sections 
shown in Fig. [31 In image (i), the black arrows point out two spots where vortices enter the film more easily. Below (c)-(j) are 
reconstructions from the magnetic data: components Jx and Jy of the current density flowing in the sample, and components 
Eix and Eiy of the inductive portion of the electric field. Color scales for x and y components are the same. Black streamlines 
of J and Ei overlay their x components. The complete set of frames is presented as the movie B JE . avi in the supplemental 
material. 
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the sample plane, 1.2 /im, the sensor tilt, 3° about the x 
axis, the strip's critical current, 1.13 A, the Hall coeffi- 
cient, 0.10 Vl/G, and the amplitude, 7.2 Oe rms, of an ap- 
plied field proportional to the applied current, explained 
below. The Hall coefficient's fit value is consistent with 
the calibrated value of 0.11 fi/G. The y position and tilt 
of the strip about the z axis are also allowed to vary, 
and the tilt is zeroed by rotating the image. These pa- 
rameters are constrained to be constant over time; the 
only change in the calculation from frame to frame is the 
(known) applied current. 

Adding the small, uniform applied field that varies 
with the applied current improves the fit, and is sug- 
gested by the sample geometry in Fig. [2fa), where the 
current returns to the right and below the bridge. Pos- 
itive returning current generates a negative field at the 
section of bridge imaged, which accords with the sign of 
the field added to the calculation. Furthermore, the ver- 
tical segment of the return lead should add a negative 
dB/dx, which is not accounted for by the calculation, 
and which does explain why the disparity between the 
data and the calculation is largest toward the lower right 
corner of the boundary. 



FIG. 3: (Color online) Cross sections of the data in Fig. (2] 
The sections are taken along the y axis at the location of the 
dashed yellow line in Fig. [2jc). The vertical dashed green 
lines mark the boundary between data and background fit. 
Successive frames in the current cycle are overlaid; their colors 
and labels match the frame labels in Fig. [J] 

Figure [2Uc)-(j) are selected from a larger set of frames, 
and represent one half-cycle of ac applied current. Ini- 
tially, the maximum current 7=1.07 A is applied and flux 
has penetrated, somewhat inhomogeneously, into both 
edges of the strip. / is apparently below the critical cur- 
rent, /c, as a central fiux-free region separates the two 
fiux fronts. This condition is necessary for our calcula- 
tion of electric field, as explained in Section IVl 

By (f), the applied current is reduced to zero, and vor- 
tices remain trapped in the edge regions. Flux of the 
opposite sign starts to enter at the edges. We then see 
this opposite flux erase and replace the trapped flux as 
negative current is applied. The succeeding half of the 
cycle approximately repeats the S, J, and E configura- 
tions shown, but with opposite signs. 

The boundary of the data is shown as a dashed green 
box in (c). The surrounding background is calculated 
with a critical state model for a thin superconducting 
strip, The main purpose of this background is to bet- 
ter match the boundaries of the image for the Fourier 
transformations described in Section IIV Al The good 
agreement between the model and the data also shows 
that much of the strip's response can be ascribed to crit- 
ical state behavior, though deviations, such as spatial 
inhomogeneity in pinning strength, are evident. 

The model is fit to the entire set of frames at once. 
The free parameters are the height of the sensor above 



IV. CURRENT RECONSTRUCTION 

A. Magnetic inversion with regularization 

The Biot-Savart law describes the magnetic field B 
generated by a current distribution J, 

B(r) ^ f / .^r'ii^^liii^. (2) 
A-K J |r-r'| 

Several authors have tackled the problem of inverting this 
relation to obtain a planar current distribution J(.T,y) 
from a planar magnetic measurement Bz(x,y)-'^''^''^ We 
tried two of these existing methods: regularization^ 
and conjugate-gradient.'^ The conjugate-gradient method 
produces current distributions dominated by unphysical 
artifacts. Further testing with simulated data suggests 
that this occurs when the current distribution extends 
outside the image boundaries, as in our images, which 
encompass only a section of the superconducting strip. 

We met with greater success using the regulariza- 
tion method with generahzed cross-validation (GCV), de- 
scribed in detail in Ref. and summarized here. Taking 
the z component of Eq. [2] and Fourier transforming in x 
and y, we find an algebraic relation between the Fourier 
transformed quantities B^^ Jx, and Jy: 

Bz{kx, ky, z) ~ —^e — {kyJx kxJy^ ■ (3) 

We assume that Jz, current flowing perpendicular to 
the plane of the film, is insignificant. We also ignore 
any z dependence of the in-plane components of current, 
approximately solving for the current density averaged 
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over the film thickness. These are reasonable approxi- 
mations in our film's geometry; its thickness, 180 nm, 
is smaller than its efi'ective London penetration depth, 
Aoff = Aafc/tanh(d/2Aafc) = 400 nm,^iS^ The specific ker- 
nel that we use to relate and Jy to Bz is that of a film 
of 180 nm thickness 

We now reduce Jx and Jy to a single unknown by not- 
ing that J must be nearly divergence-free at the operating 
frequency of 400 Hz (the resonant frequency, I/VlC, of 
a piece of strip like that imaged would be ~ 100 GHz). 
This allows us to derive both J^, and Jy from the local 
magnetization g[x,y)^ where 

J(a;, y) = -zx Vg(x, y) ^ J^ ^ dyg, Jy = (4) 

In Fourier space, 

g{kxi ky) — — (jtyJx ^ kxJy^ ■ (5) 
Eqs. [3] and [5] yield 



Thus our basic procedure is to Fourier transform B^, 
solve for 5, Fourier transform back to the real space 
g(x,y), and use Eq. |4]to obtain Jx and Jy. 

The first difficulty arises from the factor e*^^ in 
Eq.[6l Spurious high spatial frequencies (with wavelength 
greater than z, the measurement height) in the magnetic 
data are exponentially amplified in the inverted current. 
While scanning Hall probe microscopy enjoys lower noise 
and smaller z than magneto-optical imagingf2^ noise in 
our images can still dominate the reconstructed J. The 
method of regularization compensates by suppressing 
high frequencies (smoothing), and GCV determines an 
optimal amount of regularization from the data itself. In 
practice, we used GCV as a guide to choose a regulariza- 
tion parameter (A = 100 as described in Ref. 0) that we 
held fixed across the set of frames. 

Similarly, differentiating g to find the components of J 
amplifies high-frequency noise, so we use Savitsky-Golay 
smoothing to extract the derivatives. The smoothing is 
quadratic with a frame size of 2.5 /zm (5 pixels). 

The second difficulty is that the Fourier transform of 
B assumes periodicity in the vertical and horizontal di- 
rections, and mismatches between the left and right, and 
top and bottom boundaries of the B image lead to arti- 
facts dominating the image once it has been transformed, 
manipulated, and transformed back. A common solution 
is to window the data, bringing its boundaries smoothly 
to zero. This discards a large portion of each image, 
however. 

Instead, we rotate the original data (2.3° about the z 
axis) so that the bridge runs horizontally, then center it 
on a larger area (a square 256 /im on a side) in order 
to move the edge effects away from the data. We fill in 
the background with a calculated B as described in Sec- 
tion mil Only after calculating J and E^ (described in 



Section IV Ap do we crop the images back to the original 
dimensions (plus the margin seen in Fig. [2I^c)). Finally, 
before Fourier transforming, we mirror the images top- 
to-bottom in order to better match the top and bottom 
edges without going to the much larger area necessary to 
allow the field to die off, which would be more computa- 
tionally cumbersome. 

While these preparations remove artifacts associated 
with edge mismatch, we do observe a spurious bump in 
the reconstructed current at the boundary between data 
and fit, seen clearly in the cross section of Jx in Fig. [31 
However, this artifact is about ten times smaller than our 
signal, and appears to be confined to the boundary. 

B. Discussion of current density images 

The results largely agree with our expectations for a 
superconducting strip. The reconstructed current flows 
within the strip approximately in the x direction. At the 
maximum applied current (Fig.[21[c)), the current density 

i^J — \J "^x + Jy^ the edge regions of fiux penetration 
should equal the critical current density, according to the 
critical state model. J then dips down (but is not zero) 
in the vortex-free central region, as expected from the 
demagnetization effect of the strip geometryJ^i^i This 
separation into edge and central regions is clearest in 
the cross section of Jx, Fig. ^c). J averaged over the 
edge regions is 12 MA/cm^, which accords with macro- 
scopic transport measurements of Jc on a similar film (9 
MA/cm^ at 44 K).'2S We expect transport, which is sen- 
sitive to the weakest point in a superconductor, to yield 
a lower value than the spatially averaged Jc. 

The cross section of Jx also shows more current on the 
—y side, indicative of the applied field (described in Sec- 
tion lIII[) modifying the symmetric distribution one would 
expect for an applied current. While the cross sections 
vary along the length of the strip, this asymmetry is typ- 
ical. 

Jy, plotted on the same color scale as Jx, is smaller, but 
highlights where the current reroutes around apparent 
weak spots in the film. Spots producing the largest Jy 
are marked with arrows in Fig. Wi}) ■ These features in Jy 
correspond to bumps in the streamlines overlaid on Jx 
where the current spreads around these defects. In the 
following Sections, we show that these are also spots of 
high electric field and dissipation. 

V. ELECTRIC FIELD RECONSTRUCTION 

We relate electric field to the magnetic field we measure 
via Faraday's law, 

V X E = -dtB. (7) 
Taking the z component, 

dxEy ~ dyEx — —dtBz. (8) 
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This only defines E up to the gradient of a scalar. We 
therefore use the Helmholtz decomposition to separate 
the electric field into a divergence-free inductive portion, 
Ej, and a curl- free electrostatic portion, Ep (following 
the notation of Ref. |17) , 

E = E, + Ep. (9) 

Ep = — V</) where is a scalar potential and V • Ep = 
p/eo, the charge density. Our measurements determine 
Eli through Eq. [H but do not determine Ep without fur- 
ther constraints, as described in Section IV CI 

To illustrate: If we apply a dc current / > /c to our 
strip, flux flow or other resistive behavior generates an 
Ep, but we would see little time variation of the magnetic 
field when averaging over length scales greater than the 
intervortex spacing. In contrast, for I < Ic, the voltage 
and electric field (both E^ and Ep) are zero in the steady 
state. We therefore remain below /c of our superconduct- 
ing strip in order to minimize unmcasurable portions of 
Ep. 

Furthermore, below 1^, even in a dynamic state, Ep 
remains zero for a strip that is uniform in a; with no 
Hall effect. In this case, symmetry dictates that Ey — 0, 
= 0, and d,^E.j, = 0, thus V • E = 0. Then Ep is 
uniform, and zero below Ic- Thus all of the behavior we 
expect from a model strip will be contained in E^. 



To approximate the time derivative of we compute 
twice the final number of frames (200 frames spaced by 
12.5 ^s, but each still averages a 25 /zs interval) and take 
the differences between successive even frames. We can 
then calculate simultaneous electric field and current at 
the times of the odd frames. 

Once we have solved for h and transformed back to 
real space, we use Savitsky-Golay smoothing, as with J, 
to extract the partial derivatives corresponding to Eix 
and Eiy. 

Finally, Eq.jSjonly defines Eix and Eiy up to constants. 
Setting the constants is equivalent to finding a field-free 
point, to which Norris devoted much careJ^ We set the 
zeros based on the edges of the uncropped images, as 
far from inhomogeneities in the data as possible. For 
the Eix zero we use the mean of the two pixels at the 
vertical center of the strip (where we expect no E, as 
Norris pointed out) on the left and right edges. We set 
the Eiy zero to the mean of the four corners. This mean 
is zero for an ideal strip, and minimally affected by fields 
originating toward the center of the image. 



B. Discussion of Ei images 



A. Inductive electric field Ei 

To reconstruct the divergence-free E^, we proceed as 
for J, solving for a potential function h where 

Ei(x,2/) = -z X Vh{x,y) ^ E,x = dyh, E,y = -dxh 

(10) 

in which case Eq. |H] becomes a Poisson equation for h: 

~dlh~dlh = -dtB, (11) 
Vloh = dtB, (12) 

We are interested in the electric field in the sample 
plane (z = 0) rather than the measurement plane. We 
therefore need Bz{z — 0), which we obtain from the 
current distribution via the Biot-Savart law. at any 
height is easily computed by rearranging Eq. |B| to obtain 
Bz in terms of 5, the magnetization function from which 
J is derived: 

B. = ^e-'^^k-g. (13) 

It is also easiest to solve the Poisson equation in Fourier 
space, where Eg. IT^ becomes 

{ikfh^dtB,{z^Q) (14) 

and, combining Eqs. IT3l andlT4|, 



As the applied current decreases over the half cycle 
shown, Eix remains approximately zero (white) in a cen- 
tral region of the strip. Outside of this region it becomes 
negative (blue), continuing past the strip edges. This is 
also visible in the cross sections of Eix in Fig. O This 
behavior accords with our expectations for vortices mov- 
ing into the edges of the film. The central region shrinks, 
tracking the flux front as vortices enter. Although vor- 
tices and current from the previous half cycle are present 
inside the central region, the vortices remain pinned and 
therefore do not generate an electric field. 

The central region is not completely field-free, however. 
Interestingly, the field it does display — about 10 times 
smaller than the edge fields — is maximal and opposite to 
the current when dl /dt = 0, a point where our critical 
state model would dictate E — Q. Such a field would 
arise, though, from a relaxation of Jc, i.e. flux creep in 
which the vortices continue to move into the strip even 
as the current momentarily stops ramping. It leads to a 
(temporary) negative power input to the film, discussed 
in Section FVll 

Such movement while dl /dt = is also visible in the 
full set of magnetic images. While the magnetic field 
evidence alone is subject to errors in phase relative to 
the applied current, the electric field confirms that the 
relaxation is real. Such relaxation is the focus of Ref. [l3- 

Finally, the cross sections of Eix reveal an unexpected 
negative tilt, dE/dy, both in the central and outer re- 
gions, which is an error that arises from the tilt of the 
Hall sensor, as shown in Section IVHII 
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C. Electrostatic electric field Ep 

We established in Section |V] that for a uniform super- 
conducting strip with no Hall effect below /c, Ej, = 
E = Ei, even with an ac current. Indeed, the E^ that 
we observe contains all the features we expect from a 
uniform strip, as discussed in Section fVB I Our images, 
however, also reveal inhomogeneity, which could produce 
a non-zero Ep. 

Here we show that in spite of this inhomogeneity, we 
expect the total electric field E to remain approximately 
parallel to J locally, which proves sufficient to reconstruct 
Ep as the field that compensates for any component of Ej 
perpendicular to J. Our method is inspired by Ref. |l7| . 
but we do not make the additional and incorrect assump- 
tion that the component of Ep parallel to J is zerOf22. We 
also describe a rather different route to E and Ep in Sec- 
tion |VlII3 



1. Validity of the constraint E||J 

First we must justify that in our experiment, E is par- 
allel to J, emphasizing that this will not be true for all 
materials. For example, a material can have an intrinsic 
Hall effect. However, macroscopic transport measure- 
ments indicate that for YBCO in the superconducting 
state, the Hall effect is insignificant; the component of 
E perpendicular to J is at most 1000 times smaller than 
the parallel component. 

Another violation of E||J could arise from a feature 
such as a grain boundary, whose orientation prevents vor- 
tices from moving perpendicular to J. More generally, 
any gradient in superfluid condensate energy density will 
exert a force on vortices. We show here, however, that at 
least on lengths scales greater than the image resolution 
of 1 /xm, we can rule out the presence of gradients strong 
enough to compete with the pinning forces that occur on 
the scale of the coherence length, <^ « 2 nm. We put 
an upper bound on such a gradient in our material by 
assuming that the gradients in pinning strength that we 
observe stem entirely from changes in condensate energy 
density. We take J (the magnitude of the current den- 
sity) at maximum applied current as a map of Jc near 
the edges of the strip (as discussed in Section HV Bp . The 
magnitude of the pinning force per length for a single 
vortex is 

Fp = $0 Jc (16) 

from which we estimate the depth of the vortex pinning 
potential as = Fp£_. The line energy of the vortex will 
be ~ • 10, from the ratio of the calculated depairing 
current to the measured Jc-^ Thus the gradient of Fp 
yields a gradient of line energy, i.e. a force Fgrad, whose 
magnitude we compare to Fp. We find that Fp/Fg^ad 
averages 2000 over the edge regions of the strip with a 
minimum of 200. 



This analysis suggests that in our material, the Lorentz 
force from the current at Jc, along with pinning forces, 
dominate other forces felt by vortices, so that vortices 
move perpendicular to J and generate E parallel to J. 
However, we do not claim to rule out every possible ma- 
terials effect, e.g. effects that average out below our res- 
olution of 1 fim. 

2. Method of reconstructing Ep 



Starting with E^ and J, reconstructed in Sections IV Al 
and ITV| and the constraint that E = Ei -I- Ep is par- 
allel to J, we can solve for Ep. The component of Ep 
perpendicular to J must cancel that of E^ : 

Epi_ = -E,^ (17) 

where Ei±_ is calculated by subtracting from E.^ its pro- 
jection onto J. Having obtained -Epj^, we construct Ep\\, 
the component of Ep parallel to J, to satisfy 

V X Ep = (18) 

which follows from the definition of Ep. The z component 
of Eq. [18] tells us that at each point, 

d^Ep\\=d\\Ep^ (19) 

or more explicitly, 

- '^d..Ep\\ + ^dyEp\^ = ^d^Ep^ + ^dyEp^, (20) 

which is a linear, first-order partial differential equation 
for Epn{x,y). 

We tried two approaches to solving this equation for 
E'pII. In the first, we start with a trial solution, Epi, 
composed of the known Ep± and i?p|| = 0. We calculate 
its curl (V X Ep)z, which will be zero for the true Ep. 
From the curl, we reconstruct a divergence-free field F 
following the same procedure we used to reconstruct Ej 
from Bz , described in Section IV Al The reconstruction 
preserves V x F = V x Epi , so G = Epi — F is curl- free, 
as desired. However, G does not preserve Ep±. So we 
construct a new trial Epi consisting of the components 
Gil and Ep±, and iterate. The true Ep would remain 
unchanged by such a procedure. In practice, however, we 
found that the procedure misconverged to an Epi with 
larger curl than the initial trial. 

We therefore attempted a more direct integration of 
Eq.IH 

Ep\\ ^ j d\\Ep^ds (21) 

where the integral is taken along a path s that remains 
perpendicular to J at each point along its length, and 
starts at some (x,y) such that Ep\\{x,y) = 0. 

The result, constructed from many such paths, is 
shown in Fig. U) The initial data (J and E^) is taken 
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FIG. 4: (Color online) Reconstructing the electrostatic portion of electric field, Ep: All images are cropped to the area of the 
strip. The electric fields are on the same color scale, shown in (a), and the curls are on the scale shown in (c). (g)-(i) are x 
cross sections from the images above them, with locations indicated by the dashed lines of matching color. The cross sections 
also indicate the spatial scale of the images. The component of Ep perpendicular to the current, (a)-(b), is obtained from Ei 
under the assumption E||J. However, (c) shows that this component alone does not satisfy V x Ep = 0. We calculate and add 
a parallel component, yielding the complete Ep in (d)-(e). Our _Ep|| calculation requires an integration constant, defined by 
assuming Ep\\ = along the lines shown in (d) and (e). Adding i?p|| suppresses, but does not perfectly cancel, V x Ep, as seen 
in (f) compared to (c), and from the cross sections in (i). 



from frame h in Fig. ^ then cropped to the area of the 
strip. For starting points, we set Ep^^ (x, y) — along the 
horizontal line shown in Fig.Ufd). In practice, to ensure 
coverage of every pixel, we start a path at each pixel and 
work back to the zero line. For efficiency, we skip pixels 
that have been covered by previous paths. 

Ep±, shown in Fig. |lta)-(b), represents the starting 
data from which wc calculate d\\Ep±. Integration gives 
us Ep\\ , which we add (as vector components) to Epj_ 
to obtain (d)-(e). The method is far from perfect, as 
evinced by the non-zero V x Ep in (f). However, in com- 
parison to (c) , the curl is suppressed at all but the highest 
points. This is clear in the cross sections through (c) and 
(f) shown in (i). So the result of the procedure, (d)-(e), 
is closer to, but still short of, the true, curl- free Ep. 

Finally, we note that this method gives us no informa- 
tion about Ep outside the sample, where J — 0. However, 
outside the strip the charge density p/eo = V • Ep — 0. 
Then Ep — —Vcj) where (/> obeys Laplace's equation, 
— 0, with a Neumann boundary condition given 
by Ep in the strip. 



D. Discussion of Ep images 

Figure [5Uc)-(j) shows the total electric field, E = 
Ei -I- Ep, for the set of frames from Fig. [21 The com- 
plete set of frames is assembled into the movie EP . avi in 
the supplemental materiali^ Our sample clearly deviates 
from a uniform strip, in which Ep = 0. Ep is compara- 
ble in magnitude to E^, but much more inhomogeneous, 
contributing most at the spots identified by arrows in 
Fig. [2li) as weak points of the superconductor. 

Thus we demonstrate reconstruction of the total elec- 
tric field from our time-resolved magnetic images. We 
reiterate that this analysis is restricted to materials in 
which E is parallel to J, and in which we can identify a 
field-free kernel. 



VI. RECONSTRUCTED POWER INPUT 

Armed with J and E, we calculate P = J • E, the local 
power input to the film, resolved in time and space. The 
results are shown in the third row of Fig. [Sj The largest 
features are the positive edge regions where vortices move 
in as the current sweeps. We note that the instantaneous 
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FIG. 5: (Color online) (c)-(j) The total electric field, E, and power, P, for the same set of frames shown in Fig. [J] The 
images are cropped to the area of the strip, (k) Integrating P over the area of each image and the thickness of the film, 
and normalizing by the length of strip imaged, we obtain the total power input as a function of time over a cycle of applied 
current, (m) Integrating P over time, we see the spatial distribution of energy input over a cycle. Any reactive component of 
P integrates to zero, leaving only the dissipated energy. 



power shown arises from both dissipation and reactance. 
These are not easily separated (e.g. by the relative phase 
of J and i?) because of the nonlinear relationship between 
J and E. 

We can relate our local measurements to macroscopic 
transport measurements by integrating over space, shown 
in Fig. [DJk) . Each frame of the 50- frame set covering the 
central portion of the current cycle is summed and nor- 
malized by the length of strip imaged to obtain power 
per unit length, then plotted at its time within the cycle 
of applied current. These points are repeated in the first 
and fourth quarters of the cycle (in which J and E re- 
peat with opposite signs) . We note that when integrating 
over the entire sample (which we only do imperfectly by 
integrating over the image area), we expect no contribu- 
tion to the power from the electrostatic field, Ep, which, 
exerting a conservative force, cannot do work. Indeed, 
dropping Ep from the calculation shown in (k) does not 
change the result significantly. 

The power input rises as the magnitude of current in- 
creases, but then falls back to become negative as the 
current reaches its peak. An inductive response with 
zero resistance would be zero at the peak. Instead, as 
discussed in Section IV B[ this negative contribution (E 
opposite J) arises from relaxation of Jc. As the applied 
current decreases back to zero, we do not recover much 
power, as we would in a dissipation-free inductor, be- 
cause the vortex movement is irreversible — the vortices 
remain pinned. 

When integrating over time, any inductive contribu- 



tions to the instantaneous power input cancel, leav- 
ing the sum of dissipation over one cycle, shown in 
Fig. [Hl^m). As in the instantaneous power, the edge con- 
tributions dominate. Integrating this image over space, 
or equivalently integrating the power in (k) over time, 
we obtain the energy dissipation per length of conduc- 
tor, 6.6 • 10~^° J/cm/cycle. For comparison, a calcu- 
lation from the critical state model using the applied 
/ (1.07 A peak) and fit value of /c (1.13 A) yields 
1.2 • 10~^ J/cm/cycleJ^ Factors in the lower measured 
value may include suppression of high spatial frequen- 
cies in the reconstructed quantities due to regularization 
(see Section llV A|) and cropping of the image, which may 
exclude some pixels near the edges of the strip. 

This energy dissipation translates to an average power 
of 2.7 • 10"^ W/cm, which, given the applied sample cur- 
rent of 0.75 Arms, implies a voltage of 6.9 nVrms between 
the voltage taps, which are spaced by 260 /im. We have 
not attempted to verify the presence of this small but 
perhaps measurable voltage. 



VII. LOCAL E-J RELATIONS 

Over the cycle of applied current, each location in the 
film experiences a range of current densities and electric 
field strengths. By plotting E against J for each pixel, we 
obtain a local characterization of the material. For ex- 
ample, the E-J curve has often been used to test models 
of the flux pinning mechanismi^i^i^i^^i^ 
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FIG. 6: (Color online) E, the magnitude of the projection of 
the total electric field onto J, is plotted against the magnitude 
of J. Values for one column of pixels are overlaid. The color 
of a set of points indicates its position on the bridge within 
the column of colored pixels in the inset, which is an image of 
the magnitude of J at maximum applied current. The values 
for one pixel are highlighted in black and the pixel location is 
marked in the inset. 



Figure [H shows E-J curves for one column of pixels 
at the locations indicated in the inset. One curve is 
highlighted, illustrating a plausible E-J relation for a 
superconducting film. While the points at low J dis- 
play several fiY/cm of scatter, there is a clear upturn 
at Jc w 10 MA/cm^. This value is consistent with the 
12 MA/cm^ we observe in the edge regions of the strip 
at maximum applied current, and with 9 MA/cm^ from 
transport measurements of a similar film at 44 K, as dis- 
cussed in Section flVBI ^^ However, as in Section IVTl we 
caution that the electric field we measure arises from 
both dissipation and reactance, and in comparing our 
data with a purely dissipative dc E~J curve, we ignore 
reactance. 

Although many of the curves seem reasonable, they 
deviate significantly between pixels, with upturns at cur- 
rent densities ranging from 14 to 2 MA/cm^. These de- 
viations are mapped out in Fig. [7] for various electric field 
criteria, Ec- The value of each pixel is the lowest J for 
which Ec is exceeded. If Ec is not exceeded, the pixel 
is plotted as white. If we use Fig. Elja), which shows 
J at maximum applied current, as an estimate of Jc, 
it seems that many of the curves hit Ec at erroneously 
low J. However, the spatial variation at least partially 
reflects genuine inhomogeneity of the material, because 
the points of low Jc along the edges of (b)-(d) correlate 
with low points in (a). 
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FIG. 7: (a) The magnitude of J at maximum applied current 
compared to Jc extracted from E-J curves for various electric 
field criteria Ec, shown in (b)-(d). All are on the same color 
scale, shown in (b). Pixels that did not reach Ec are white. 
Though (b)-(d) show more scatter, areas of low Jc generally 
match those of low J in (a). 



VIII. RECONSTRUCTION ERRORS 

Throughout the paper we have tried to point out physi- 
cal assumptions and where they may break down. Even if 
our assumptions hold, however, the reconstruction proce- 
dure can introduce errors and amplify uncertainties. The 
magnitudes of these effects are difficult to predict ana- 
lytically because of the many complex numerical trans- 
formations involved. Instead, we estimate the errors by 
executing the procedure on simulated data, for which we 
know exact solutions to compare to the reconstructed 
quantities. Figure [5] compares cross sections from the 
resultant exact and reconstructed images. 

We generate the data with the same critical state 
model used for the background of the magnetic images, 
as described in Section |TTT1^-"^ The parameters are the 
same, with the applied current (-0.02 A) corresponding 
to Fig. [2If). We choose this frame because it includes 
regions of zero current, which present a worst case for 
reconstructing Ep, as described below. 

As the input to the reconstruction procedure, we cal- 
culate an image of the magnetic field in the measurement 
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FIG. 8: (Color online) Errors introduced by the reconstruction procedure are analyzed by executing the procedure on model 
data. The exact quantities, known from the model, are plotted as blue lines. From them, the magnetic field in the measurement 
plane and its time derivative are calculated, noise is added, and the quantities are reconstructed. The red dots are y cross 
sections through the centers of the reconstructed images. The residuals are the reconstructed values minus the exact values. 
The only variation as a function of x comes from the noise, so the exact values of the y components are zero. 



plane (z=1.2 ^m). This is identical to the background of 
Fig. [HJf), but extends over the entire image area. As 
in the background, we account for the sensor's tilt (3° 
about the x axis), which adds a small portion of —By to 
the measured "-Bz." We then add normally distributed 
noise within the measurement region (the subset of the 
image defined by the dashed green box in Fig.[2Jc)). The 
amplitude of the noise matches that observed in the real 
data outside the strip. 

Using the reconstruction procedure, we obtain J^, Jy, 
and Bz{z = 0), shown in the first column of Fig. [H The 
reconstructed Jx smooths the sharp corners of the true 
Jx- It also displays 0.03 MA/cm^ rms of noise, as does 
the reconstructed Jy . The only variation as a function of 



X comes from the noise added to B; the exact values of 
all y components (and Epx) are zero. 

To reconstruct E, we similarly calculate B for ap- 
plied currents corresponding to 12.5 /xs before and after 
Fig. lUf). From each B image we reconstruct Bz{z — 0), 
and subtract to approximate dBz{z = 0)/dt, shown in 
Fig. IHI^k) . We then follow the reconstruction procedure 
for Ei (as in Section lyX)) and Ep (Section [VC]). 

We see in Fig. [SJe) that an erroneous negative slope is 
present in the residual of the reconstructed Eix . This er- 
ror can be traced to the small portion of —By in B, which 
is not accounted for by the reconstruction procedure. In 
principle, the procedure could be adapted to assume a 
specified linear combination of By and B^, though re- 
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lations such as Eq. [6] would become more complicated. 
One could also calculate an approximate By from the re- 
constructed Jx, subtract it from the measured B, then 
iterate, reconstructing a more accurate and Bz each 
iteration. This procedure has been successfully applied to 
removing in-plane field components from magneto-optic 
images. '^^ Both procedures require precise knowledge of 
the sensor tilt. 

Although small compared to the error in E^^ , other er- 
rors are worth noting. For one, we compare the y compo- 
nents of the inductive electric field and total electric field 
(Fig. [Sfh) versus (i)). The addition of Epy suppresses 
Ey (from about 6 • 10~^ to 2 • 10~^ /zV/cm) everywhere 
except near y = ±18 /im, where Ey fluctuates with x 
by 0.1 /iV/cm rms over the image. At these locations, 
J is approximately zero, making the direction of J com- 
pletely uncertain, which renders Eqs. [TT^El inaccurate. 
Furthermore, such errors may scale up with the larger Ep 
present in the real data. Thus we must admit ~ 100% 
uncertainty in Ep near regions in which J drops below its 
noise level, 0.03 MA/cm^. More generally, we point out 
that this error analysis is performed for a homogeneous 
model, which may not account for some features of real, 
inhomogeneous systems. 



Alternate reconstruction of E from Ei 



where 

d J — d Jx 



^dtBz 



{dxJ)Jy - {dyJ)Jx 

(27) 

are known. The solution to this differential equation is2& 



P{J) 



J^' ui J)q{ J)dJ 
u{J) 



(28) 



where u = exp ^ p( J)dJ^ . We use the boundary con- 
dition p(0) = 0. 

So the ingredients are the reconstructed dtBz{z = 0) 
and J, along with its spatial derivatives. These quantities 
are arranged by time; we re-sort them by J. Then the 
integrals with respect to J are calculated by the trape- 
zoid method. Applying Eq. [^5] at each pixel, we obtain 
p{J) and therefore E(J) at each point. We could then 
reconstruct Ep = E — if desired. 

In practice, our results are dominated by noise. In 
our images, the spatial derivatives of J are of order 
10"'^ MA/cm^ while their uncertainties (based on the 
0.03 MA/cm^ uncertainty in J) are several orders of mag- 
nitude larger. This method may be effective, however, for 
other samples or measurement techniques. 



Where there are large uncertainties in Ep, one can use 
an alternative route to reconstructing E that bypasses 
Ep. This method starts with V x E, to which Ep (which 
is curl free by definition) does not contribute. As in Sec- 
tion |VCl we assume E is parallel to J, in which case we 
can write 

E = p(a;,y,J)J (22) 

and describe E via the scalar p (the local resistivity). 
Swapping sides and taking the curl, we have 

(Vp) X J + p(V X J) = V X E (23) 

Taking the z component and applying Faraday's law 
(Eq. [8]) to the right-hand side, 

{dxP)Jy - {dyp)Jx + p{dxJy - dyJx) = -dfBz (24) 

To proceed, we assume that the spatial dependence of p 
arises solely from its dependence on J, i.e. that the E-J 
relation is constant over the single-pixel scale at which 
we solve this equation. Then, 

^ [idxJ)Jy - idyJ)Jx]+pidxJy-dyJx) = -dtBz (25) 
which we can write in a standard form 

^+pp{J)^q{J) (26) 



IX. CONCLUSION 

We have shown that time-resolved magnetic imaging of 
a superconducting thin film yields a complete characteri- 
zation of its electromagnetic properties, including distri- 
butions of current flow, electric field, power, and local E- 
J relations. We also point out the physical assumptions 
and requirements behind the mathematical transforma- 
tions. The technique is compatible with various methods 
of magnetic imaging. It requires time resolution commen- 
surate with the sample conditions being studied, but this 
capability is not limited to scanning Hall probe micro- 
scopes. Magneto-optics can use high-speed frame grab- 
bing or a phase-locked short pulse technique to acquire 
similar data sets (albeit in a different pixel sequence) 
Alternatively, the average response can be acquired from 
each pixel in succession, as is done here, at the expense 
of longer total acquisition time. 

An applied field, Bz{t), can be substituted for our ap- 
plied current, and its time dependence need not be si- 
nusoidal. For example, the analysis is applicable to the 
transient response following a change in applied current 
or field, as in Ref. [l3. However, by demonstrating the 
technique on a superconducting strip carrying an ap- 
plied current at 400 Hz, a realistic operating frequency 
for power applicationsji^ we point out that an important 
use of this work is the characterization of material inho- 
mogeneity and its effect on ac losses. Such imaging may 
also provide an efficient method for measuring the effects 
of complex sample geometries/^ 
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As mentioned in Section IVIIi E-J curves give insight 
into the flux pinning mcchanism^i^i^i^^i^ By resolv- 
ing the variation in E versus J with position and time, 
this analysis may allow one to dissect the behavior of a 
heterogeneous sample, correlating pinning dynamics with 
materials properties. One may also examine correlations 
with quantities, such as magnetic field strength, that vary 
with space or time, fully accounting for the heterogeneous 
self-field experienced by different parts of the sample. 
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